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Abstract 

In this paper we explore the structure and applicability of the Dis- 
tributed Measurement Calculus (DMC), an assembly language for dis- 
tributed measurement-based quantum computations. We describe the for- 
mal language's syntax and semantics, both operational and denotational, 
and state several properties that are crucial to the practical usability of 
our language, such as equivalence of our semantics, as well as composition- 
ality and context-freeness of DMC programs. We show how to put these 
properties to use by constructing a composite program that implements 
distributed controlled operations, in the knowledge that the semantics of 
this program does not change under the various composition operations. 
Our formal model is the basis of a quantum virtual machine construction 
for distributed quantum computations, which we elaborate upon in the 
latter part of this work. This virtual machine embodies the formal se- 
mantics of DMC such that programming execution no longer needs to be 
analysed by hand. Far from a literal translation, it requires a substantial 
concretisation of the formal model at the level of data structures, naming 
conventions and abstraction mechanisms. At the same time we provide 
automatisation techniques for program specification where possible to ob- 
tain an expressive and user-friendly programming environment. 



Keywords: quantum computing, measurement-based quantum computing, 
distributed computing, formal models, virtual machines. 

1 Introduction 

During the last decennia, quantum information has managed to become a sig- 
nificant field of research in the exact and applied sciences. Although it is a 
relatively new discipline one can currently discern several sub-disciplines such 
as quantum cryptography, information theory, computability, error correction. 
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fault tolerance, computations and of course there is also the experimental re - 
search on the construction of quantum computers ( Nielsen and Chuand . 2000( ). 
Nevertheless, the development of quantum information as a proper computa- 
tional domain of computer science is lagging behind. Indeed there is no such 
thing as a quantum computational paradigm. By this we mean a framework in 
which quantum problems can be expressed and solved in terms of data struc- 
tures, algorithms, techniques such as abstraction, all of these wrapped up in 
an associated programming language. Paradigm building has proved to be an 
extremely useful approach in computer science. It has led to theoretically equal 
but practically very different programming frameworks, such as functional, im- 
perative, logic and object-oriented programming. For this reason we expect this 
approach to be crucial also in developing quantum programming paradigms. 

The first step in paradigm building is to construct a low-level quantum 
programming language and determine its properties. By low-level we mean that 
we need to define syntactical expressions for each physically allowed quantum 
operation: preparation, unitary transformation, measurement, combined with 
classical control expressions if so desired. The syntax in its own is not the goal, 
but rather a means by which to facilitate investigations with computer science 
techniques. First, we need to determine the functionality of a quantum program, 
its semantics. The most obvious way to do this is the operational semantics, the 
most practical is the denotational semantics. While in the former a program's 
meaning is given as a sequence of state-changing operations, in the latter it is 
instead a mathematical object. Paramount is linking both, so that one can use 
whichever in future analyses. Through a programming language's semantics 
one can investigate notions such as composition and context-freeness. These 
properties are crucial when one wants to build more complex programs. Indeed, 
they allow the semantics of these larger programs to be built up from that 
of smaller components using rules for composition, rather than having to be 
determined from scratch. While these properties may seem obvious, computer 
science is littered with examples where they were mistakenly taken to be true, 
leading to problems in prog ramming language development (see for example 
( Brock and Ackermanl . 198lh ). 

Recent advances in quantum communication and cryptographic computa- 
tions motivate the need for a programming paradigm centred on distributed 
quantum computations. In a distributed system one has concurrently acting 
agents in separated locations, operating locally on a quantum state, which may 
be entangled over agents, and coordinating their actions via communication. 
Formal frameworks for distributed quantum computation have only very re- 
cently begun to appear. Typically, these are a combination of classical pro- 
cess theory, which formalises notions of concurrency and communication, the 
quantum circuit model, i.e. local operations are unitary transformations of an 
age nt's qubits, and given ini t ial shared entanglement. First, there is the work 
in ( Lalire and Jorrandi 2004 : Jorrand and Lalire . 2005 ). which is directly built 
upon classical process calculi. While this model profits from being closely related 
to existing classical models, the disadvantage is that it is hard to focus on quan- 
tum behaviour. A different approach is advanced in the model known as com- 



muni cating quantum processes or CQP ()Gav and NaearaianL l2004t iGav et al. 



20051), where the typical communication primitive s of proc e ss cal culi are com- 



bined with computational primitives from QPL (jSelingeii 120031) . The basic 



model of CQP is more transparent. This model serves as a basis for the dcvcl- 
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opment of formal verification techniques, in particular for proving the security 
of larger scale quantum communication systems. In related work, a probabilistic 
model-checki ng tool bui l t upo n an existing automated veri fication component 
is de veloped ( Gav et all . 12005 ). There is also the work in ( Adao and Mateus . 



20051), which is specifically tailored to security issues in cryptographic proto- 



cols. In our work, however, we are much more interested in the expressiveness 
of quantum distributed computations and the behavioural properties of dis- 
tributed primitives. In a way, we take the inverse approach, assuming that 
computations are well-defined and investigating what programming concepts 
are at work, instead of the other way around. 

While the fact that formal verification tools for distributed quantum compu- 
tation are currently under development may suggest that a mature distributed 
paradigm already exists, this is actually not the case. Distributed protocols 
are still very much conceived on intuition, and there is no good notion or for- 
malisation of the programming concepts that are at work. We therefore require 
adequate formal tools with which to explore and evolve the distributed quantum 
computing paradigm. In this paper we explore the structure and applicability of 



the D istributed Measurement Calculus (DMC) (jD'Hondtl . 120051: iDanos et al. 



2005[ ). an assembly language for distributed measurement-based quantum com- 



putations. We describe its syntax and semantics, both operational and denota- 
tional, and state several properties that are crucial to the practical usability of 
our language, such as equivalence of our semantics, as well as compositionality 
and context-freeness of DMC programs. We show how to put these properties 
to use by constructing a composite program that implements distributed con- 
trolled operations, and demonstrate that the semantics of this program does 
not change under the various composition operations. Finally, we elaborate on 
an implementation for the DMC language which we developed under the form 
of a virtual machine, a platform-independent programming environment that 
abstracts away details of the underlying hardware or operating system. This 
virtual framework is a crucial step in providing DMC as an experimentation 
platform for distributed quantum computations, as reasoning within the formal 
model by hand proves quite cumbersome for ev en small coniputat ions. At the 



basis of our work lies the measurement calculus ( Danos et a^.l . 120071 ) . a low- level 



quantum programming language for measurement-based quantum computations 
from which we explore the distributed paradigm. Next to the obvious advan- 
tage of starting from a proper formal framework, measurement-based models are 
well-suited as a starting point for distributed quantum computations because 
they are inherently distributed. What is important here is that the well-known 
physical framework of quantum computation is ported to an equivalent com- 
puter science framework, opening up the field towards investigations from this 
branch of science as well. 

The structure of this paper is as follows. In the next section we discuss the 
syntax of our language, while Sec. [3] covers the semantics of DMC and lists its 
properties. Sec. S] applies the previous to a concrete example, a composite proto- 
col implementing distributed controlled gates. We discuss the implementation 
of the quantum virtual machine for DMC in Sec. [SI and conclude in Sec. [HI 
Some basic knowledge of quantum computation is assumed; for the rea.der no t 
familiar with the domain we refer to the excellent ( Nielsen and ChuanBi[2000l) . 
Our approach in this article is to explain the notions of the model by example, 
rather than providing a series of formal definitions which are hard to interpret 
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and for which space it too hmited. However, we stress that the model is a rig- 
orous one, and, while this paper is a stand-alone document, refer the interested 
reader to t he appendix for the full formal semantics of the DMC language and 



(jP'Hondtl . l2005t iDanos et all 120051) for complete definitions 



2 Syntax 

The language we propose is, broadly speaking, an assembly language for dis- 
tributed measurement-based quantum computations. It is an assembly language 
in that it provides syntax only for the most basic operations while at the same 
time being universal. It is measure ment-based as we bu ild our language on top 



of the measurement calculus (MC) (jPanos et aZ.I . l2007f ). an assembly language 



for measurement-based quantum computations. The latter depart from the 
usual circuit-based approach to quantum computing, where unitary transfor- 
mations are the driving force of computations. While measurement operations 
were long seen as a necessary but disruptive part of quantum computing, in al- 
gorithms such as teleportation they act as an essential part of the computation. 
This gave rise to models where the computation is steered by pre-established 
generic entangl ement combined w i th me asurements, such as the one-way quan- 



tum computer (jRaussendorf et aZ.1 . 120031 ). Because measurements are inherently 



probabilistic, correction operations are required that are conditionally applied 
depending on previous measurement outcomes, thus rendering the computation 
deterministic. All of this is nicely captured in the measurement calculus. For 
this reason, as well as the inherent distributive aspect of measurement-based 
models, MC is an ideal basis from which to develop our language. 

We describe our language model in three layers. The base layer consists of 
MC patterns^ which describe local quantum computations. These are combined 
with distribution primitives into the notion of agents in the middle layer. Fi- 
nally, agents and their shared entanglement resources are bundled into networks 
in the top layer. 



Patterns. In MC a pattern is defined by a sequence of commands together 
with sets of qubits for working, input and output memory. As an example 
consider the following pattern, which, as we explain below, implements the 
Hadamard operation. 



7^(l,2):=({l,2},{l},{2},X|Wf£;i2) • (1) 

All qubits arc uniquely identified using numbers. The first argument denotes 
that the pattern has a computation space of two qubits, referenced by 1 and 
2. The next two arguments specify the pattern's inputs and outputs. Working 
qubits that arc not input qubits arc initialised to |+) = |0) + The last 
argument is the pattern's command sequence, a list of operations taken from 
a basic set and executed from right to left, analogous to matrix applications. 
Subscripts denote qubit arguments of the operation, while corrections arc con- 
ditionally executed depending on their superscript. Concretely, one runs the 
pattern H by preparing the first qubit in some input state \ip) and the second 
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in state |+), then entanghng both qubits with the controhed-Z operatiorQ to 
obtain CZi2i\'4') ® |+)) (syntax: E12). Next, the first qubit is measured in the 
{!+), |— )} basis (syntax: M^), where |— } = |0) — |1). Finally, an X correction is 
applied on the output qubit if the measurement outcome was si (syntax: X^^). 
Here si defines a signal — simply or 1 - coming from the X-basis measurement 
on qubit 1 in . If the signal is 0, the operation is not performed. Measure- 
ments are considered to be destructive, which is why qubit 1 does not appear 
in the output set. 

A general pattern is denoted 7-'(V, /, O, A), with computation space V, inputs 
/ and outputs O (together called the pattern type), and command sequence A 
that consists of entanglements Eij, measurements *[M°'Y, or corrections Xf or 
Zf, where «, j e 1/, a G [0, 27r] and S Z2. Allowed measurements are those 
in the XF-plane of the Bloch sphere, and are specified by the angle a {M^ is 
actually syntactic sugar for M^). Measurement angles may also be conditioned 
by signals, written *[Af"]'*, with (— 1)*q; + tir being the actual measurement 
angle. The four basic inst ructions together with signal conditioning su ffice to 



make the model universal (iRaussendorf et all l2003t iDanos et a/.l . 12007 ) 



Patterns can be combined into larger ones to create arbitrary quantum com- 
putations. The sequential composition of patterns Vi ~ (Vi, /i, Oi, ^1) and 
'P2 = (V2,/2, 02,^2), with Oi = h, is defined as 

V2V1 := {VilJV2,h,02,A2Ai) , (2) 

while the parallel composition of the same patterns is defined as 

Vi®V2 ■■= (Fl Wl^2,/l W/2,Ol WO2, A^2) • (3) 

Note that one can always rename qubits for parallel composition to become 
possible, while sequential composition also needs I2 and Oi to have the same 
cardinality. 

As an example, here is the pattern to create a 3-fold GHZ-state |000) + 

GHZu3 = ({1, 2, 2, 3, 3}, ■, {1, 2, 3}, H(3, i)E^^n{2, 2)E,^) , (4) 

where 2 and 3 are working qubits. 

MC is equipped with a powerful standardisation theorem which provides a 
procedure for bringing any pattern, such as the one above, into EMC-iovm, 
i.e. with entanglements first and corrections last. This is important from an 
experimental implementation perspective, and also corresponds nicely with the 
typical structure of a distributed quantum protocol where one starts out with 
a shared entanglement resource. In fact within MC we can already express 
distributed computations such as teleportation, 

({l,2,3},{l},{3},X3^^Z3^iM2^Mfi?i2i?23) . (5) 

While the pattern describes the intended computation, we find no notion of 
separate parties participating in the teleportation, and neither of the required 
communication between them. Purely by convention, and because everybody 
knows the protocol, we can say that qubits 1 and 2 belong to Alice and qubit 



controUed-Z operation on two qubits applies the Z operation to the second qubit 
provided the first is set to 1. 
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3 to Bob. It is still hard to see if and what both parties have to communicate, 
e.g. Bob needs si and S2 but has no access to qubits 1 and 2, thus Alice needs 
to perform the measurements and communicate the results to Bob. We want 
this information to be explicit and obvious in the language, which in turn makes 
it easier to describe and reason about distributed quantum programs. To see 
that this is important, try to identify what is implemented with the following 
distributed protocol, 

({0,0,l,l,2,2},.,{0,l,2},Z2^^Zi^iZo^oM5^jf^i?ooi?iii?22) • (6) 

Here M^^^ is a pattern for GHZ-measurement. Note that pattern (0]) imple- 
ments the unitary transformation between the diagonal basis and the GHZ-basis 
+ Hence a GHZ-measurement is executed by applying the inverse 

pattern followed by a measurement in the diagonal basis. We will get back to 
this pattern in Sec. 01 

Agents. An agent embodies a single computation node running in isolation 
in a distributed algorithm, i.e. a processor or a party such as Alice or Bob. 
Each agent has a local command sequence, which operates on a set of resources 
contained within its environment. Agent instructions arc either measurement 
patterns or communication primitives. For example, the teleportation pattern 
([5]) really needs the following agent definitions for Alice and Bob, respectively, 

A: {l,2}.d.S2.c\si.M^M^Ei2 and B : {3}.Xl^Z^\c7s2.c7si . (7) 

As we can see an agent definition consists of a type in curly braces, which spec- 
ifies what qubits an agent owns, and an instruction sequence. Next to local 
quantum operations, written in MC language, we specify the exchange of the 
signals si and S2 between Alice and Bob via the communication primitives c? 
and c!. What is missing here as compared to pattern (O is the prior shared 
entanglement E23, as indeed it is impossible to factorise this part of the pro- 
tocol into agent definitions. There is no other option than to specify shared 
entanglement separately in the network definition for the full protocol, as we 
shall see below. Of course the qubits are still local to the agents - as reflected 
in the type - only their description is not. 

Formally, an agent is defined by an expression A(i, o) : Q.£, where the type 
Q is a set of qubit references and £ is a finite instruction sequence composed 
of pattern command sequences A, classical message reception els and sending 
c!s, where s is a signal, and quantum reception qc?g and sending qclq, where 
g is a qubit reference. Corresponding communication actions are synchronised, 
meaning that an agent executing a receive will pause its program until the 
required agent has sent a message on the same channel and vice versa. Classical 
inputs i and outputs o are used in protocols such as superdense coding, and are 
not mentioned when there are none. Also, working qubits required by patterns 
are initialised to |+) as before and are not specified in the type Q. 

Networks. A network of agents consists of several agents that execute their 
command sequence concurrently. Typically, quantum resources are shared be- 
tween agents - indeed most distributed quantum protocols benefit from some 
type of shared entanglement being present prior to the start of the protocol. We 
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have no way of splitting the representation of these states over agent definitions, 
and instead have to specify agents separately in a network's definition. This is 
why a network specification is more than just a collection of agent definitions, a 
feature specific to quantum concurrent frameworks and absent in classical con- 
currency. To keep track of how a shared resource is distributed each agent's 
type Q specifies which qubits of a shared resource are local to that agent. 

At this point we can finally write down the full tclcportation protocol, in- 
cluding all relevant distributed aspects. 

TP := A : {1, 2}.c\s2.c\si.M^ E12 

B : {3}.X;^^Z^i.c?S2.c?si (8) 

II £^23 

Analogous to process algebraic notation we use a vertical bar [ to separate 
concurrently executing agents. Shared network resources are specified after the 
double bar ||. For teleportation this is the state £^23, which was produced and 
handed out to Alice and Bob prior to the execution of the protocol. If so desired 
the establishment of shared resources can itself be seen as a distributed protocol 
involving a server agent polled for entanglement services. We now likewise lift 
the pattern given in to a distributed setting, obtaining the following 

ES := L : {0, 1, 2}.co!so.ci!si.C2!s2-A^oif ^ ■ 
Ao : {0}.Zo"«.co?xg 

Ai : {l}.Zl\c^lxi (9) 
A2 : {2}.Zl\c2lx2 

II -^00-^11-^22 > 

where M^^^ is a pattern executing a GHZ-r aeasurement. This is the the entan- 



glem ent swapping protocol for three agents (jZukowski et all Il993t iBose et al. 



1998I ). which produces a GHZ-state shared between Aq, Ai and A2. Because 
the resulting GHZ state is in the diagonal basis we have Z rather than X cor- 
rections in the network specification. 

Arbitrary networks are denoted A/" = \iAi{ii,Oi) : Qi.Si \\ a. For simplicity, 
we assume that networks of agents satisfy a number of definiteness conditions 
which ensure that the computation is well-defined. For example, an agent may 
only operate on qubits he owns, and every communication event needs to have 
a corresponding dual event in the network. Of course these issues are important 
but we are glossing over them here as we assume checking these occurs in a 
pre-compilation step rather than at runtime. 

Our extensions to the measurement calculus have made distributed notions 
explicit. We can now see directly from a protocol specification what agents 
have to communicate to whom as well as which assumptions are made about 
non-local entanglement resources. Also, because each agent's instructions are 
expressed separately, we do not impose any particular execution order for local 
quantum operations. This makes it clear which parts of the protocols can and 
cannot run concurrently. 
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3 Semantics 



In the previous section we established the syntax of DMC. The next step is 
to estabhsh the meaning of a program written down in this language, i.e. its 
semantics. Formal semantics is a means of assigning objects to chunks of code 
so that one can reason with these objects rather than with the code itself. In 
this section we describe an operational as well as a denotational semantics for 
DMC programs and, with this in hand, state some important properties of our 
framework. 



Definition. A network's operational semantics reflects how it affects the state 
of a distributed system on which the network is run. That is, we associate with 
each agent a local state (quantum and classical), and specify how the network 
specification updates these. Since we always have to take the shared entangle- 
ment resources into account the operational semantics does not decompose into 
state transformers for each of the agents separately. On top of this a crucial in- 
gredient of the semantics specifies how quantum resources move throughout the 
system. Concurrency is not really an issue for the semantics because, due to the 
linearity of quantum mechanics, the order in which local operations are applied 
is unimportant. Hence we pick an order arbitrarily and derive the semantics for 
this case. 

Concretely, the operational semantics of a network is found by collapsing 
individual small-step semantical rules into one transition system. The full small- 
step semantics of DMC can be found in the appendix; it consists of a number of 
quantum mechanical rules - corresponding to measurement pattern commands 
- and a number of communication rules. These affect the global quantum state 
{EMC commands), the local resources (quantum communication), and the local 
memory {M commands and classical communication). Since measurement is 
probabilistic, so are the small-step rule s. The formal semantics of measure ment 
pattern commands is well established ( Danos et all . l2007t D'Hondt , 2005 ) , and 



just needs to be lifted to the DMC setting. We elaborate here the semantics of 
typically distributed concepts introduced in the previous section. The formal 
semantics of MC relies on the component F, the outcome map, which contains 
a number of bindings from signal names to measurement outcomes. In DMC 
the outcome map is lifted to a local agent memory recording also the values 
of classical messages. Concretely, a classical communication event between two 
agents has the following semantics, given that Ft, the local memory of Bob, 
evaluates the name y to the value u, 

A : {Qa, ra).Sa-c7x \ B : (g,, F,).£:,.c!y 
=^A:{Qa,Ta[x^v]).£a\B:{Qb,rb).£b ■ 

Here Ta[x ]. v] means that we assign a new binding of a; to w in the local 
memory of Alice. These rules are typically much harder to read and write 
down formally than to apply concretely. What we are saying here is that if a 
classical communication takes place, the value is looked up by the sending agent, 
and received and bound to a new name by the receiving agent, after which 
the computation {£a and Eh) continues. Sending and receiving qubits over a 
quantum channel changes the types Q of the agents involved. An agent sending 
a qubit can no longer perform any operations on it, therefore the corresponding 
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qubit reference is removed from while it is added to the receiving agent. 
Formally, the rule for the exchange of a qubit with reference i is given by 



A : ((5a,ra).fa-qc?« | B I (Qfe , Tf,) fc .qc !i 

=^A:(Q,U{z},ra).fa |B:(Qb\{i},rb).£fc . ^ ' 

It goes without saying that the way in which the actual communication of a qubit 
is implemented is more intricate than this simple rule. However, by providing 
quantum communication as primitive syntax we are precisely choosing not to 
get into these implementation matters. Both rules mentioned here do not affect 
the global quantum state, which is why it is not mentioned. This is not the 
case for the small-step semantics for pattern commands, which moreover may 
be probabilistic. The full small- s tep se mantics can be found in the appendix or 
in (iD'Hondtl . l2005t iDanos et ad l2005h . 

We obtain the operational semantics by defining computation paths in the 
usual way and gathering all small steps in a computation path in one big step 
from initial to final state in that path. As such the operational semantics of 
a quantum distributed network is essentially a probabilistic transition system 
(PTS). However, since resource allocation is crucial, we have to augment this 
PTS with information on how qubits move throughout the network. This is 
formalised as a type signature. We denote the operational semantics of a net- 
work M by |A/']op. For example, the operational semantics of TP is given by 
the deterministic transition system 

[rPlop:({l},-)-^(-,{3}).pi=>P3 , (12) 

where p is the density matrix specifying the input quantum state to be teleported 
and subscripts indicate qubit systems. The operational semantics of the ES 
network given in Q is also deterministic. 

lESU ■■ ■) ^ {0}, {!}, {2}) . =^ GHZ^^^ , (13) 

where the superscript D denotes that the resulting state is in the diagonal basis. 
Note that we only write real quantum inputs in the type while not mentioning 
entanglement resources. The good thing about operational semantics is that in 
principle it can be derived automatically by induction on the small-step rules. 

Denotational semantics is a second means of assigning a formal meaning to 
a chunk of code, this time by way of mathematical objects. If we look at the 
skeleton of a distributed protocol, i.e. the equivalent non-distributed pattern, 
we find a multi-local probabilistic quantum operation, which is mathematically 
represented by (a special type of) completely positive map (CPM). Since proto- 
cols with different distribution of resources are observationally different, we have 
to pair this CPM with a function mapping input to output resources, formally 
represented by a type signature. We denote this type of semantics by lA/Jde for 
a network M . For example, the denotational semantics of TP is given by the 
map 

[TPl,e:({l},-)^(-,{3}).I, (14) 

i.e. TP implements the identity map from qubit 1 to qubit 3. Here we see the 
importance of the type signature in specifying the semantics. The subtle differ- 
ence between both types of semantics becomes more apparent once one starts 
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investigating more complex protocols involving mixed states and probabilistic 
semantics. 

As we see from the examples both types of semantics are quite similar, and 
indeed we have the following result, which we state without proof. 

Proposition 3.1 There is a precise correspondence between the operational 
and the denotational semantics of networks of agents, that is to say 

VM, A/'2 : M =op J\f2 ^ Ml =de M2 ■ (15) 



Two networks are semantically equivalent if they have the same semantics!! 
While the equivalence between semantics may seem obvious (certainly from the 
example) , it is crucial to prove this statement formally, as it allows one to switch 
between semantics where appropriate without giving the issue further thought. 
We note that in general this equivalence is not guaranteed. 

Building larger protocols. With syntax and semantics in hand we can now 
consider constructing larger networks from smaller components. We consider 
parallel as well as sequential composition of networks, denoted by ® and o 
respectively. These operations are defined in the obvious way so we will not spell 
out concrete definitions here. Suffice to say that one needs to pay attention to 
agent as well as qubit names to ensure that networks are connected as desired. 
More importantly we need to make sure that these forms of composition are 
consistent with the semantics, as is indeed the case. 

Proposition 3.2 The semantics of networks is compositional, i.e. 

[AA2oA/-il-IAA,Io|A/-iI . (16) 
|A/l0A/-2] = IA/-2|®IA/-il (17) 



Here we have a first application of Prop. 13.11 as first we do not need to specify 
which type of semantics we mean in the statement of Prop. 13.21 , and second we 
can choose the most convenient semantics for the proof, which is the denota- 
tional semantics. While the compositions on the left hand side are at the level 
of agent and network definitions, the composition of network semantics on the 
right hand side is a functional one at the level of typ e signatures, PTS's and 



CPMs. For the full proof we refer to (jD'Hondtl . 120051 ): an example is given in 
Sec. HI 

A second important notion when constructing larger protocols, closely re- 
lated to compositionality, is that of contexts. Informally speaking a context is a 
program with a "hole" in which a network specification can be inserted, typically 
denoted C[-]. On some occasions, one finds that programs that arc considered 
equivalent in isolation no longer behave in the same way when placed within the 
context of a larger program. This is particularly the case in concurrent systems . 



A historical example is the Brock- Ackerman anomaly ( Brock and AckermanL 



^For example, one can prove that qubit communication between two agents is semantically 
equivalent to the teleportation network. 
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by which it was reahsed that simple input-output behaviour is not enough 
to reason about equivalence of concurrent systems. In our framework contexts 
C[A/'] are given by arbitrary network compositions C in which our network M 
is placed. What is different with ordinary compositionality is that the quan- 
tum resources of the network A/" are no longer considered to be independent of 
those of its context. Indeed in Prop. quantum resources are considered to be 
provided independently, such that composite networks operate on disentangled 
inputs. While this is sensible when each of the networks operate in isolation, it 
is less so when a network is only one factor in a complex compositional struc- 
ture. Hence we first need to show that our semantics is independent of so-called 
entanglement contexts. This is a consequence of the following proposition, which 
we state here without proof. 

Proposition 3.3 Suppose C : pA — > Pb = X^fe LkPAL\ ^ completely pos- 
itive map. Then for all quantum states pAc applying C to the A-part of pAc 
results in 

PBC ^^{Lk®Ic)pAc{Ll®Ic)- (18) 

k 

The proof, though easy, is not trivial. Using this proposition and composition- 
ality we have the following important result. 

Corollary 3.4 Equivalence holds in arbitrary contexts, that is, if Mi = A/2, 
then for all network contexts C[-\ we have that C[A/i] = C[A/2]. 



4 Applications 



We now show how all of the formal ingredients can be put to use in a concrete 
example. The distributed primitive we will investigate is that of a distributed 
remotely controlled gate ( Yimsiriwattana and Lomonacol 2005 ). Near- future 
quantum computers are expected to have only a limited number of qubits per 
machine. Even in quantuni simul ating environments the current qubit limit is 
only about 36 (jRaedt et aU 120061 ). Hence one can imagine that quantum com- 
putations need to proceed much like cluster computations today, with resources 
spread over different processors in order to make them feasible. One common 
situation would be where the central processor needs to execute a controlled 
operatio n with the target qubits spread over a group of agents. Following the 
ideas in ( Yimsiriwattana and Lomonacol 20051) , once we have a GHZ resource 
we can execute a distributed controlled gate. Calling the central processor L (for 
Leader), and assuming there are two subordinate processors A and B, the trick 
is to establish a shared control qubit between target agents, which is achieved 
through the share control (SC) protocol as follows. 



SC 



L : {c,0}.ci,2!so. Mq^-Boc 
Ai : {l}.Xf°.ci?xo 
A2 : {2}.X2"".C2?xo 

I GHZ^^2 



(19) 



Here qubit c is the input control qubit which is in the state a|0) -I- /3|1). That 
the protocol indeed establishes its goal may be seen from its semantics, which 
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can be derived unambiguously to be 

[5Cl„p:({c},.,.)-^({c},{l},{2}) 

(a|0)+/3|l))c=^(a|000)+/?|lll))ci2 . ^ 

Once this type of shared entanglement is in place each target agent just has to 
execute a local controlled unitary gate with as control its qubit in the shared 
entanglement resource. Note that because we have context-independence of the 
semantics, target qubits may be entangled over different agents. However, in 
order for the distribution approach to be possible at all, the controlled unitary 
CU must have U = Ui ® U2 where, Ui and U2 operate on a number of target 
qubits smaller or equal than the maximum available qubits for each agent. 

We see that the SC protocol requires GHZ entanglement. It is probably 
realistic to assume that in a quantum network each agent can ask for Bell-state 
entanglement with the central server. It is not so realistic to assume that groups 
of agents can demand direct GHZ entanglement whenever they need it; rather, 
we expect GHZ entanglement to be produced via the entanglement swapping 
protocol. Note that, since we need GHZ entanglement between L, Ai and A2, 
we need to compose L with Ao in the ES protocol as presented earlier in 
that is 

L := AooL : {0, 0, 1, 2}.Zo"'\ci!si.C2!s2.Mo°if'' • (21) 

What we are actually doing here is composing the SC protocol with no shared 
resource with the ES protocol to establish the resource, and our semantics 
ensures that this is something we can do unambiguously. Indeed, we have 

ISCoESj^lSCjofESj 

= ({c},.,.)^({c},{l},{2}) (22) 
. (a|0) + m)c =^ (a|000) + /3|lll))ci2, 

with |5C] and {ESj given in ^ and ^ respectively 

More agents. The networks that we defined for implementing the distributed 
remote gate protocol can be generalised to n agents. This requires generalised 
procedures for GHZ-measurement, entanglement swapping and establishing a 
shared control qubit. First, an n-fold GHZ state is produced by generalising 
the 3-GHZ-pattern given in Q. That is, through the pattern with no input 
qubits, output qubits {1,2,..., n} and an event sequence interleaving E and H 
operations, as follows, 

Gi/Zi...„ = n{n, n)E^,,_,)f, . . . E^^n{i, S)E^^n{2, 2)E,^ , (23) 

where the hatted qubits are again working qubits. Again, a GHZ-measurement 
is executed by applying the inverse pattern followed by a diagonal-basis mea- 
surement. This leads to the pattern 

M^3f = M^M? . . . M9E,.^n{2, 2)E^^H{3, 3)E^;, . . . E^,,_,)fMn, n) , (24) 

with input qubits the qubits {1, 2, . . . , n} to be measured and with no outcome 
qubits. Using this sub-pattern we can establish GHZ-entanglement between the 
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leader L and agents Ai through to An by the genera lised entanglement swap- 
ping protocol ( Zukowski et at . 1993t Bose et ai . 1998 ). which has the following 



network specification, 
ES := 



L:{OA.-.,n}.Z^\ic,\sj)-^,.M^"J' 
\U A. : {^}.Zf^Ci?x, 
II ®^oEu- 



(25) 



This is just the generalisation of network ^ where we have merged agents L 
and Aq for the purpose of establishing a shared control qubit as before. The 
signal s = Sq . . . Sfi corresponds to a projection on the GHZ-statc \s) + \ s). This 
network has the following semantics, 



[iS^lop :(•,..., •)^({0},{1},...,W).0 



GHZ, 



D 

0...n 



(26) 



After establishing GHZ-cntanglemcnt between agents a shared control qubit is 
obtained through the following protocol. 



SC:^ L:{c,0}.Cilsa.M^Eoc ■ 
\U A, : .cjxo 
II GHZq ^ 

The semantics of this network is given by 

I5Clop:({c}, )->({c}, {!},..., M) 

(a|0)+/3|l))e=^(a|0)®("+^^+/3|l)^ 



(27) 



■(n+l) 



)cl 



(28) 



Control qubit c can now indirectly control unitary operations at the sites of 
all agents by having agent A,; execute a local pattern for a controlled unitary 
gate where the control is its qubit i and the targets are locally available qubits. 
Again, for this approach to be possible the controlled unitary CU must have 
U = Ui ® . . . ®Un where each of the sub-unitary Ui is executable by agent A^ . 



5 Virtualisation 

In the above we introduced a formal language for distributed quantum compu- 
tation. Providing a number of tools for constructing higher-level programs, it 
should be seen as a first step in a bottom-up construction of a distributed quan- 
tum programming paradigm. However, the previous section clearly shows how 
cumbersome it is to describe and evaluate computations purely within the formal 
model. Rather, these developments are only really useful when one thinks of the 
language DMC in terms of a quantum virtual machine (QVM). A virtual ma- 
chine is a platform-independent programming environment that abstracts away 
details of the underlying hardware or operating system. In our setting it is a 
low-level language abstraction layer which executes DMC programs independent 
of the actual implementation of quantum operations, which could be executed 
by any of several existing quantum simulators or even by a physical quantum 
computer. As a mediator between a set of low-level basic quantum gates and 
the construction of more complex quantum programs, a QVM forms a crucial 
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layer in a tiered quantum computation architecture ( Svore et ai . 2006[) . It is 



especially valuable if one is concerned with developing higher-level distributed 
quantum computation languages through experimentation and abstraction, for 
which a user-friendly and flexible programming environment is invaluable. 

Translating the formal model into a virtual machine has the obvious bene- 
fits of automating program execution and composition. Furthermore, low-level 
inspection during a pre-compilation step can automatically determine a number 
of issues such as well-defincdncss of code. While the formal language we have 
discussed in the above provides the backbone for our virtual machine, its actual 
implementation is far from trivial. Issues such as well-definedness of programs, 
naming of variables within local computations as well as in larger program con- 
texts and efficiency all come into play in a more concrete sense that is absent 
in the more abstract formal model. In what follows below we give details on 
our QVM implementation. We first discuss the QVM for sequential, MC-based 
computations, then extending the platform further towards a distributed ver- 
sion thereof. We use the adjectives formal and virtual to differentiate between 
similar objects in the formal model and the virtual one whenever the context is 
unclear. 

5.1 The quantum virtual machine 

The first step in the development of an execution environment for the mea- 
surement calculus is translating formal measurement patterns into structured 
data for the c omputer to work with. We represent this data with symbolic or 
s-expressions 



(lMcCarthvl . [l960l) . as popularised by the programming language 



Lisp. At the same time we chose to adhere insofar as possible to the nota- 
tion used in the formal model. Before discussing the full virtual syntax below, 
we give the example of the Hadamard pattern from Eq.(IT|), which is expressed 
virtually as 

n-.^ ?o) (?i) (?o) ((E ?o ?i) (M ?o 0) (X ?o (s ?i)))) . (29) 

Similarly to the formal setting we see a pattern expressed as a list of three qubit 
sets {V, I, O) and a command sequence. Question marks in qubit names indi- 
cate that they are variables, subject to renaming and instantiation to concrete 
qubit references. Important to note is that in the command sequence operations 
are executed from left to right, in the formal model notation this was the re- 
verse. This change is due to both implementation reasons and computer science 
tradition. 

Our current execution environment is split in two separate layers: an exe- 
cution and composition layer. The execution layer takes a command sequence 
as an input and is in charge of performing the quantum operations it specifies. 
This command sequence is obtained after the data representation of a mea- 
surement pattern is assembled into the low-level language that the execution 
layer can understand - essentially by replacing qubit variables with concrete 
qubit references. Hence in terms of abstractness this layer lies below the one in 
which ordinary patterns are defined. The composition layer, on the other hand, 
adds an abstraction layer on top of pattern definition in terms of the pattern 
composition structures from MC. It compiles the composition of patterns into 
a single pattern data representation, by merging arbitrary compositions of mul- 
tiple patterns into a single new pattern. In the rest of this section we discuss 
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implementation details of each of these components, focusing on respectively 
the execution layer and its pattern assembler, as well as the composition layer 
and its compiler. 



5.1.1 Execution layer 

One arrives at the execution layer whenever a pattern has to run, either directly 
as a result of its definition being called or indirectly after a composition of 
patterns was compiled to an expression ready for execution. The execution 
layer consists of an interpreter which directly executes each operation in its input 
language. This interpreter is the core of our quantum virtual machine (QVM), so 
called to stress the low-level and machine-generated nature of its input language. 
It is essentially an automatic version of the operational semantics described 
earlier in Sec. |3] and specified in full in the appendix. The input language 
syntax for the QVM is obtained by turning MC notation into an s-cxpression 
form, specified as follows in a BNF notation: 



<sequence> 

<instruction> 

<correction> 



<entanglement> 
<signal> 



( { <instruction> } ) 

<correction> I <measurement> I <entanglement> 
= ( X <quref > [ <signal> ] ) I 
( Y <quref > [ <signal> ] ) 
<measurement> ::= ( M <quref> <angle> [ <signal> ] [ <signal> ] ) 
= ( E <quref > <quref > ) 
I 1 I <input-name> I 
( s <quref> ) I ( + <sigiial> { signal } ) 

(30) 

where the symbols in boldface indicate syntax (round brackets and literals), 
while square and curly brackets are meta-syntax symbols which denote op- 
tion (zero or one occurrence) and repetition (occur any finite number of times) 
respectively. The only real difference with the formal syntax is the notation 
( s <quref > ) for signals, referring explicitly to measurement outcomes rather 
than relying on the naming convention s„ to denote the outcome on qubit n. 
The execution layer interpreter takes any command sequence written in this 
language and executes it according to the semantics of MC. Note that in this 
way we already obtain a low-level quantum programming language. For exam- 
ple we can write and execute the Hadamard gate with the following command 
sequence: 



((E 1) (MO 0) (X 1 (s 0))) . (31) 

Transforming a pattern definition Eq. ([29| to its assembled version Eq. ((3T|) is 
relatively simple. Concretely, the pattern assembler replaces each qubit vari- 
able name in a pattern's command sequence by concrete qubit references. For 
every distinct variable name used we choose a unique integer and replace every 
occurrence of the variable by it. At this stage only the input qubit set is used, 
in particular for initialisation of the (non-input) auxiliary qubits to the default 
|-|-)-state. As we shall see below both input and output set qubit sets are used 
when combining patterns in the composition layer. 

Practically, the interpreter takes besides a command sequence also an envi- 
ronment, which contains the current state of execution. After an initialisation 
step the first operation in the sequence is evaluated, after which the interpreter 
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recursively calls itself with an environment modified by this operation. Just as 
in the formal model, the environment for the interpreter consists of a quantum 
as well as a classical part, respectively denoted Q and F below. They are defined 
as follows, where T stands for tangle. 



r 

Q 

T 



= {o,i) (32) 
= T\T®Q (33) 

= ({'7l,'?2, • • -grj, \^qiq2...q„)),n G N (34) 



Differences with the formal model arise from practical considerations. The clas- 
sical state r is split into the outcome map o and the input map i to simplify 
lookup of qubit and input names respectively. Effectively, evaluating an expres- 
sion of the form (s n), where 71 is a number, becomes a lookup in o. In the 
formal model the quantum state (called p in the semantics) is presented as a 
single qubit state by creating a tensor product of all computation qubits during 
initialisation. Given the computational explosion of operations on larger ten- 
sors, this is highly impractical in a classical simulation environment. For this 
reason we introduce the concept of a tangle, a set of qubits which is known to be 
disentangled from the full quantum state and thus allowing a more compact rep- 
resentation. Hence in the QVM the full qubit state Q is composed of a number 
of individual tangle objects T , which are, evidently, updated during execution. 
Concretely, during initialisation a new tangle Ti : {{qi\, \'ip)) is created for each 
qubit, where a single qubit reference, and auxiliary qubits have = \+). 
After interaction due to an entanglement operation Eq.q. qubit references qi and 
qj are put into the same tangle T^... : ({gi, g^, . . . }, \'4'q^qj...)), while the original 
tangles Ti and Tj are destroyed. Here the dots represent other qubits which 
have interacted with qubits i or j earlier in the computation. After initialisa- 
tion each command in the command sequence is executed closely following the 
semantics of the formal model, albeit updated to the slight variations in auxil- 
iary structures as exhibited above. In our proof-of-concept implementation, we 
have implemented the semantics directly using numerical linear algebra. This 
implemen tation was bas e d on an Lisp-based simulation environment developed 



earlier in (jPesmet et al\ . 120061 ) 



5.1.2 Composition layer 

The language understood by our execution layer is essentially MC without com- 
position operators. Writing larger quantum algorithms like this by hand quickly 
becomes a tedious and error-prone process. The measurement calculus intro- 
duces measurement patterns and pattern composition which simplify this task. 
Pattern definition enables local reasoning while composition allows the creation 
of more complex programs out of smaller parts. In this way more general pat- 
terns can be defined by using parameters, and algorithms can be programmed by 
composing them out of smaller building blocks. All these abstractions are taken 
care of by the composition layer which we explain now. Its compiler effectively 
machine-generates command sequences in the language we see above. 

Patterns are command sequences where all qubit names are categorised as ei- 
ther input, output or working qubits. Defining patterns is a straightforward pro- 
cess: the command sequence syntax is the same as the execution layer language. 
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Figure 1: Graphical representation of AA"s pattern composition. 



but instead of concrete qubit references variable names may be used. An exam- 
ple is shown in Eq. (P^ . On top of this a pattern definition can be parametrise d, 



which is necessary for patterns such as the J'(a)-pattern (jPanos et al\ . \200Tu . 

By using qubit variables we facilitate the rewriting of concrete qubit names, 
which in turn allows patterns to be composed by concatenating their command 
sequences while sharing the correct qubit names. Composing two patterns into 
a larger one is defined in Eqs. ([5]) and These definitions rely on the pro- 
grammer to prepare the composition by renaming qubit names and, in case of 
sequential composition, tensoring identity patterns where needed to make sure 
the condition n V2 = Oi = I2 is met. In a setting where we wish the pro- 
grammer to be able to compose arbitrary patterns in intricate configurations 
it is crucial to automatise this renaming process. Our pattern compiler does 
precisely that. It allows the programmer to declare his intentions in non-trivial 
composition cases, whereas a standard rule is applied in absence thereof. We 
explain the general composition case first, as we define the standard rule in 
terms thereof later on. 

As an example let us consider the controlled- not pattern AX, which is defined 
in terms of the H and AZ = ({1, 2}, {1, 2}, {1, 2}, E12) patterns, as follows. 

AX (1(1) ® H(3, 4)) o AZ(1, 3) o (1(1) ® ^(2, 3)) (35) 

where T is the identity pattern, used as filler to match the number of in and 
output qubits. This composition pattern simply follows from the analogous 
matrix identity. The same intent can also be expressed without using concrete 
qubit names as follows 



AX {{I{qs)<E>H{qe,q7)) o AZ{q5,q4) o (^(qg) «, H(<7i , 92)) 
with q2 = 94, 93 = 95, 94 96,95 = 9i)}, 

which is just the textual version of Figure [1] (in fact, as we shall see below, 
our framework also provides a graphical tool for specifying constraints, much 
as in the figure). By simply matching variable names, we can derive Eq. (j35p 
automatically from Eq. ((36|) . In fact we can simplify Eq. (p6)) further by getting 
rid of identity patterns and only expressing non-trivial matching qubit names, 
which is specified in our (left-to-right) syntax as follows, 

{('H(9i,'72)),AZ(95,94),H(96,97)),{(92,94),(94,96)}} ■ (37) 

Note that the syntax uses tuples for constraints, and we have reverted to our 
left-to-right evaluation order as everywhere in the virtual model. 
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In our approach the programmer has to hst patterns in order and give pairs 
of qubits to specify how patterns are to be composed. Note that if no pairs 
are specified this corresponds qubits not being linked up, or in other words 
to a parallel composition rather than a sequential one. One can then compile 
pattern composition by instantiating the involved measurement patterns with 
fresh variable names and subsequently matching the correct variable pairs. Now 
that we have this generalised notation for pattern composition in place, we 
can derive an automatic renaming procedure corresponding to the standard 
rules for composing patterns (e.g. linking up the first pattern's output qubits 
one by one with the second pattern's input qubits for sequential composition). 
This allows us to automate the tedious part of the qubit renaming process, 
while at the same time reducing the chance for human error on larger pattern 
compositions. However, at the same time more arbitrary composition structures 
may still be defined through the qubit pair syntax provided to the programmer. 
Nevertheless, one should view the latter as a more advanced use of the QVM, 
while in most cases qubit bindings are generated as per the standard case. 

Once a set of qubit pairs has been specified for a certain pattern composition 
- either automatically or by the user himself - the automated renaming process 
works by generating a set of bindings that map the qubit variable names in 
every pattern in the composition to new ones, such that the chosen names are 
equal if they appear in the same pair. This is essentially the same process as one 
is assumed to execute manually when composing patterns in the formal model. 
To ensure that the construction of new bindings occurs in a well-defined and 
finite manner, elementary compositions in a composite structure are processed 
in topological order. To be precise, an arbitrary pattern composition expression 
such as the one in Eq. (j37|) . is viewed as a graph, with patterns as nodes and 
qubit variable pairs as edges. Since we allow only pattern combinations that 
form directed acyclic graphs, there is a unique topological ordering on the list 
of nodes, and this is the order in which we evaluate composition bindings. The 
full set of bindings B is constructed iteratively through the following rules on 
the pairs of every elementary pattern composition in topological order. 

i?(g) = B(g') = 0,9c fresh ^^^^ 



iq,q')^B[q,/q][q,/q'] 

Bjq) = qc 

iq,q')~^B[qJq'] 

B{q') = qc 



(39) 

{q,q')^B[qc/q] ^^^^ 

In other words, when both names in the pair do not appear in the binding list, 
rule (|55)) will trigger. A fresh qubit variable name qc is chosen and added as 
binding for both variable names in the pair. Rules ((39|) & (|40|) ensure that if a 
binding already exists for one of the variable names in the pair, the other will 
use the same binding. The topological sort will ensure that at all times only 
one of the three rules will execute. 

Once the binding set B is constructed qubit variables in each of the compos- 
ing patterns are substituted with their bound value, P/ = {B{q)\\/q £ Pi}. After 
this renaming process all patterns in the composition can finally be merged pair- 
wise by joining qubit sets in the right way and appending command sequences. 
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Since our composition mechanism encompasses both definitions in the formal 
model our rules for joining qubits sets are also slightly more complicated (as 
they are more general). Concretely, we have the following definition, where 
patterns are assumed to have passed the renaming process already. 

Definition 5.1 The composition of patterns Vi = (Fi, /i, Oi, ^i) and V2 = 

(V2, 12, O2, A2) is defined as the pattern V ~ (Vi U V2, O, A1A2) where 



Indeed, qubit variables fi:om the old input and output sets that were matched 
become auxiliary qubits and hence arc no longer represented in the final input 
and output set. 

We stress again that except in the most complex cases, the programmer does 
not have to explicitly denote qubit bindings when creating pattern compositions. 
Indeed, in most situations these can be created automatically or one can use 
implicit notation for qubit bindings as in the formal model. For example, we 
automatically derive regular sequential composition as specified in Eq. ^ when 
both patterns have distinct and unique variable names while the number of out- 
put and input qubits are the same. The automatic renaming process ensures 
that the condition I2 = Oi holds by automatically creating a composition ex- 
pression of Pi and P2 while identifying each qubit in Oi with the corresponding 
one in I2 (i.e. in sequential order). Likewise, parallel composition as specified 
in Eq. ^ can be trivially derived, simply by not specifying any qubit pairs. 
Without being renamed, the qubit names of each pattern in the parallel com- 
position remain unique. Because of this, Eqs. (|4T|) - (|42|) essentially merge the 
corresponding qubit spaces in the right way. 

With these shortcuts for the original composition operations in hand, we may 
now return to our example pattern AX. Indeed, while one way of specifying 
this composition is given in Eq. (j37p . by relying on the shortcut for sequential 
composition we can also express this composition without qubit names as: 



The compiler will instantiate each pattern with fresh variable names, which, 
since they are unique, will precisely lead to the desired compositions as specified 
in the expressions above. 

5.2 The distributed quantum virtual machine 

Now that our virtual machine for sequential MC computations is in place, we 
move on to its extension into a distributed version for the DMC language elabo- 
rated earlier in this article. Incorporating distribution into our framework comes 
down to extending each of the abstraction layers. Our philosophy has been to 
view network definitions as essentially specialised forms of pattern composition. 
That is, they are dealt with in a generalisation of the composition layer from 
Sec. 15.1.21 and compiled towards a set of distributed pattern definitions, i.e. 
patterns which also allow communication commands. This set of distributed 
patterns, together with information on how the network is set up in terms of 



/ = /lU(/2\0l) 
0= (Ol\/2)U02 . 



(41) 
(42) 



AX {I®-H)o AZo{I®'H) . 



(43) 
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shared resources and agent channels agents, is then assembled to an executable 
form and sent to a distributed extension of the execution layer from Sec. 15.11 
We discuss both the execution layer and the composition or network layer of 
the distributed quantum virtual machine (DQVM in short) in Sees. 15.2.11 and 
15.2.21 below, first giving a quick overview here. 

In the formal model agents are defined through a pattern with distribu- 
tion extensions, while the shared resources may be viewed as a regular non- 
distributed pattern. Before discussing the full virtual syntax below, we give the 
example of the leader agent L in the entanglement swapping protocol in Eq. ^ . 
Its data representation in our QVM is as follows: 

L := ((?qo, ?qi, ?q2) (?cho, ?chi, Ycha) 

(Aff^f4 V,, (send ?cho (s ?qo)) (send ?chi (s ?qi)) (send Ychj (s Tqs)))) . 

(44) 

As before, an agent pattern consists of a qubit sort, specifying which qubits it 
owns, and a command sequence in which pattern commands as well as commu- 
nication operations (here send instead of ! ) can be used. The main difference 
we see is that agents now also have a channel sort, in line with the resource 
sensitivities of distributed computations and with our naming mechanisms for 
variables, be it qubits or channels. 

A network definition is essentially a grouping of multiple agent patterns (the 
single bar |) together with a shared resource pattern (the double bar ||). In our 
setting we have to augment this with a so-called network configuration, a list of 
qubit and channel variable pairs to organise composition of the network com- 
ponents in the right way. That is, qubit pairs specify how the shared resource 
pattern needs to compose with the different agents patterns, as in Sec. 15.1.21 
while channel pairs indicate how communication between agents is organised. 
Concretely, a network definition is compiled into a list of multiple command 
sequences that are to be executed concurrently. The network configuration will 
enable the automatic renaming process to match the names of various qubit and 
channel names in these command sequences. We explain this process in more 
detail in Sec. l5.2.2l below. 

5.2.1 Execution layer 

The language understood by the execution layer of our distributed quantum vir- 
tual machine (DQVM in short) is the lowest-level language in our framework. 
It is arrived at after a list of distributed patterns, compiled from a network def- 
inition at a higher abstraction layer, is in turn assembled to a list of command 
sequences devoid of variable names, much as in Sec. 15. II Each of these sequences 
is a data representation of an agent's command sequence in a distributed net- 
work, and the interpreter for DMC needs to run these sequences concurrently. 
Note that the shared resource pattern (behind the double bars 1 1 ) is dealt with 
during a compilation step executed prior to arriving at the execution layer, so 
that qubits are initialised in the right way once we are ready for actual exe- 
cution. While the agent abstraction is not explicitly present in the execution 
layer, it does provide the required distributed functionality, namely channel 
communication operations and concurrent execution of multiple sequences. 
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We extend the input language syntax for the QVM (specified in Eq. ([50]) ) 
to support DMC computations as follows: 



<network prograiii> 
<agent sequence> 
<agent instruct ion> 



( { <agent sequence> } ) 
( { <agent instruction> } ) 
<instruction> I 
<channel-send> I <chaainel-receive> 
<channel-send> ::= ( send <channel-name> <signal> ) 

<channel-receive> ::= ( recv <channel-name> <input-name> ) . 

(45) 

This is again a mirror of the formal syntax, bar the use of the labels send and 
recv instead of ! and ? to indicate message sending and reception respectively 
- this to avoid confusion with variable names. Also communication commands 
now carry an extra variable specifying the channel to be used for communication, 
an improvement over the formal model where the use of channels was somewhat 
implicit. Finally, we see no syntax for the entanglement resource, which is 
instead viewed as part of the initialisation procedure carried out in the network 
layer of our virtual machine. That is, during compilation from the network layer 
we first run the pattern to construct the shared resources through the QVM. 
The quantum state resulting from this execution is then passed on to the initial 
environment of the execution layer. 

The main functionality of the DQVM's interpreter over the QVM's is to 
do with communication operations and the scheduling of multiple command 
sequences. Practically, this means that the current state of execution retained 
in the interpreter's environment not only contains information on the quantum 
and classical states of the network {Q and F respectively) but also on issues to 
do with channel usage and scheduling. We call the latter the network state N 
and add it to the interpreter's environment. The network state = {C, P} 
holds the channel map C and the network program P. The latter is a simple 
collection of each agent's command sequence at each point in the computation, 
and is used by he interpreter to keep track of each command sequence as it is 
dynamically scheduled for execution. Each command sequence is executed to the 
point where it is empty or a send or receive operation blocks. A next command 
sequence will then take turn via a round robin selection and be executed in the 
same fashion. This round robin system ensures that sequences with blocking 
operations become available for execution at later time, successfully executing 
the send or receive operation in question if a matching receive or send has been 
performed by a previous sequence's turn. This process halts when all command 
sequences are empty. Deadlocks or similar errors will cause the interpreter to 
get stuck in a loop. We note that the formal semantics of a network is the 
same for any chosen schedule and for this reason we may choose a schedule at 
will without h aving to worry abo ut the network's behaviour deviating from the 



intended one (|Danos et all 120051 ) 



The second component of the network state is the channel map C. It simply 
maps channel names to values (sent along the channel) and is used to implement 
the communication primitives in a straightforward way: 

C(ch) = F(s) = v 
^ ' ^ ' (46) 



F, C, ((send c/i s) f) ^ T,C[v/chli£) 
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C{ch) =v 

— . (47) 

r, C, ((recv c/i name) f) — >■ T[v /narae\,C[0 / ch],{£) 

Here ch is a channel name, T{s) denotes the value of the signal s with respect 
to the given outcome map and £ is the rest of the execution sequence. In other 
words, the send operation associates the value to be sent to the appropriate 
channel name in the channel map C, while the corresponding receive operation 
removes the value from C and stores it in the classical state. In case where 
a channel already contains a value during a send operation along a particular 
channel, the send operation blocks, and likewise for a receive operation on an 
empty channel: 

C{ch) = V 

C, ((send c/iu') £) ~ C,((send ch v') £) ^"^^^ 



C{ch) = 

{(xecv ch name) £) — >■ {(recv ch name) £) 

Note that these rules constitute a concrctisation of the semantical rule for clas- 
sical communication in the formal model (see Eq. llOp . Here we have chosen 
a concrete semantics using a form of blocking, in the sense that execution of 
a particular command sequence may halt until certain requirements are met, 
embodied in the equations above. 

Next to the network state the interpreter's environment also keeps track of 
the quantum and the classical state of the network. For practical reasons we 
reuse the original QVM's computation state as described in Sec. I5.1.1[ i.e. we 
have one global quantum state Q and one global classical state F for the entire 
network. This simplifies the structure of the DQVM, allowing it to focus on 
the network-specific features while the QVM is used for the execution of regular 
MC operations. Indeed, the quantum state Q is passed through unchanged to 
the QVM which executes regular MC (i.e. quantum) operations. The execution 
schedule described above imposes a synchronised access to this global state 
(Q, F), while our naming procedure ensures that names arc unique so that there 
are no clashes between classical variable names of different agents. Note that a 
more sophisticated implementation would adhere to the DMC's formal semantics 
where each agent has its own local state. For local quantum states this is to 
some extent captured by our use of tangles for disentangled quantum states 
(cfr. Eq. ([5^ 1. For classical states this would require a separate environment 
for each agent such that a received value is stored in the local classical state of 
the receiving agent. 



5.2.2 Network layer 

The network layer is the place where the programmer defines distributed pro- 
grams, such as for example the network for entanglement sharing specified for- 
mally in Eq. (|9]). A network consists of a list of agent definitions, a shared 
resource pattern, and a network configuration which specifies how resources are 
distributed and connected. Concretely we need to specify a list of patterns, 
namely the extended patterns of all involved agents and a single regular pattern 
for the shared resources. For example, to define the ES network from Eq. ^ we 
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need the leader agent pattern given in Eq. HJ] together with the agent pattern 
A and the resource pattern TZ: 



where V = I = D = (?a ?b ?c ?d ?e ?f). The fuU network is specified by 
the expression ES (7?., {L, A, A, A}, A^C) , where NC is the network con- 
figuration. As in Sec. 15.1.21 each of these patterns is instantiated with distinct 
qubit and channel variable names. The network configuration specified drives 
the composition of patterns in a network, which steers an automatic renaming 
process much as in Sec. 15.1.21 This network configuration contains two things. 
First, a list of qubit variable pairs linking output qubit variables of the shared 
resource pattern to input qubit variables of some agent pattern. Second, be- 
sides qubit variables there are now also channel variables to be matched between 
agent patterns, due to the communication extension to their command sequence. 
However, matching of channel names can be performed by the same renaming 
process defined earlier. 

After the automated renaming process has finished, agent patterns are as- 
sembled and collected in a list, forming the multiple command sequences list 
that is used as input to the DQVM's interpreter. The pattern for shared re- 
sources does not appear in the command sequences list for the DQVM. Instead, 
it is assembled and executed by the QVM as an the initialisation step, and 
passed on as initial quantum state Q in the environment to the DQVM when it 
starts executing the agent programs of the network. 

5.3 A graphical user interface 

In the above we gave an overview on the design and architecture of our virtual 
machine for the DMC language. While the subject of this article has to do 
with the inner structure of our framework, the goal of the latter is neverthe- 
less to provide a user-friendly programming environment. For this reason we 
have developed a graphical user interface (GUI) on top of the virtual machine 
to facilitate experimentation. This GUI tool currently supports only regular 
(non-distributed) patterns, allowing definition of patterns from scratch as well 
as in terms of compositions of existing patterns. Figure 15.31 shows the GUI 
being used to define a pattern, the W3 entanglement pattern, as a composition 
of known patterns; at the bottom we find the compiled version of that same 
pattern. As seen previously in this section our framework allows for a general, 
more explicit form of composition, such that explicit set of qubit name pairs 
subsumes the implicit method of matching qubits by manual rewriting. While 
expressing these pairs in writing is still a feat for large programs, a graphical 
notation similar to Figure [1] is much more expressive. Our GUI uses a similar 
notation where black boxes denote in- and outputs and the user may connect 
these in arbitrary ways to concretise the desired connections. This technique 
allows a quantum programmer to express complex quantum algorithms more 
easily, certainly for composition structures with non-trivial connections between 
component patterns. 



A : 



7^ : 



((?q) (?ch) ((recv ?ch (s ?v)) (X ?q (s ?v)))) 
(V I ((E ?a ?b) (E ?c ?d) (E ?e ?f ))) , 



(50) 
(51) 
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Actions 
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Saved Patterns' 
Pattern Composition Editor 




C Compile Patterns j 
Measurement Calculus Evaluator 



CM 78 C- 833 
CE 78 773 
CE 65 783 
CX 78 C5 3S33 
CM 38 C- 833 
CE 38 7B3 



Figure 2: Graphical User Interface of the QVM's design tool, showing the com- 
position and compilation of the 3-qubit W entanglement state preparation pat- 
tern. 



6 Conclusion 

In the above we describe an assembly language for distributed measurement- 
based quantum (or DMC) computations in all its aspects. While the first half 
of this article deals with the formal model, the second half elaborates on a 
virtual framework developed in close relationship with the formal model, i.e. a 
programming environment for the DMC language. 

DMC programs satisfy several formal properties crucial to the practical 
usability of the language, such as compositionality and context-freeness. We 
showed how to put these properties to use by formally implementing a compos- 
ite program to control operations in a distributed setting, and demonstrated 
that the semantics does not change under the various composition operations. 
DMC was developed with expressiveness in mind, a crucial property when estab- 
lishing the first layer in a distributed quantum programming paradigm. Indeed 
only through experimentation combined with abstraction can one hope to move 
towards a higher level in the language hierarchy. Our first experiments on pa- 
per already prove that DMC is indeed capable of expressing more complicated 
programs and that the formal features are necessary and sufficient in deter- 
mining their functionalit y. This is very different in flavour from earlier formal 
frameworks in this area ( Gav and Nagaraian , 2004 : Jorrand and Lalir^ . 2005 ). 
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which are much more concerned with issues such as verification, and focus on 
providing tools such as type systems to facilitate investigations of this nature. 

While the first half or this article shows that it is possible in principle to write 
DMC programs within the formal model, the limits of its usability are acutely 
felt even for the relatively small example of distributed controlled operations. In- 
deed, as an experimentation platform our formal framework is only really useful 
in terms of a quantum virtual machine (QVM) , a programming environment for 
the DMC language automating execution and composition of DMC programs. 
In the second half of this work we discuss a first implementation of such a virtual 
machine, a proof-of-concept execution environment effectively virtualising the 
formal DMC model. The benefits of having a QVM are many: as a first layer 
of a tiered quantum computation architecture ( Svore et all l2006l ) . a platform 
for automated verification and mo del- checking, and maybe most importantly, 
a basis from which to develop higher-order distributed quantum computation 
languages through experimentation and abstraction. Our QVM is built up in 
terms of several abstraction layers, most importantly a platform-independent 
execution layer to deal with the low-level semantics of basic patterns, and a com- 
position layer to create and compose larger programs. The composition layer 
comes with an associated compiler that translates any compositional structure 
to one single pattern definition. Any pattern definition, be it specified directly 
by the user or produced by the compiler, is assembled into an expression that the 
execution layer understands, essentially a command sequence where all variable 
names have been replaced by concrete references in the intended way. While the 
execution layer is defined independently of any actual implementation, we used 
a Lisp-based simulation environment in our experiments to evaluate execution- 
layer expressions. The QVM is extended with distributed structures at all layers 
into a distributed quantum virtual machine (DQVM) which allows specification 
and execution of arbitrary distributed networks. Finally, a graphical user inter- 
face (GUI) is added to facilitate usability of the framework. We note that due to 
the fact that some aspects of the formal model necessitate further concretisation 
in a virtual model, there are some differences in syntax and semantics between 
the two. For example the implicit naming conventions for variables, channels, 
and in the composition of programs, required a concrete design in the virtual 
setting. 

The virtual machine developed in this article, while covering almost all as- 
pects of the formal model, is a first implementation and as such there are several 
avenues for improvement. We list these here moving from higher-abstraction 
layers to lower ones. First, the distributed layer requires several extensions to 
be fully compatible with the formal model, most importantly by giving agents 
a more prominent role and allowing network composition. In order to develop 
the GUI into a fully-featured development tool, it also needs to be lifted to 
a distributed setting: adding agents and networks to the graphical notation 
is top priority. Ultimately the goal is to make the graphical notation into a 
self-contained language, allowing the programmer to specify arbitrary programs 
without needing to create ex-nihilo patterns, which involves writing command 
sequence code by hand. Second, much work can be done at the level of optimi- 
sation of the execution layer, which now relies on conventional rather than opti- 
mal data structures for its implementation. We are in the process of developing 
tailor-built data structures, so enhancing the performance of our framework by 
relying on domain-specific measurement calculus optimisations. For example. 
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we are currently looking into the stabiliser calculus as a means to efficiently exe- 
cute MC operations where possible, switching to a different representation when 
operations unsupported by the stabiliser calculus are performed. Optimisations 
at the level of command sequences are also possible, such as detecting and di- 
rectly initialising cluster states rather than constructing them incrementally by 
executing the required entanglements. On the other hand, classically simulating 
large entangled states means dealing with an exponential blow-u p of computa- 



tiona l time and space. For this reason the so-called EMC form (jDanos et al 



2007t ) (which puts entanglements first) is not always the preferred one in a sim- 



ulation setting, since we need to minimise the size of entangled states during 
the length of the computation. Finding the right balance between direct cluster 
state generation and exploiting classical resources at their fullest is an exer- 
cise which is currently underway. Finally, we are heavily looking into parallel 
computing techniques to improve the simulation of quantum operations, at the 
moment carried out by straightforward linear algebraic techniques. Concretely, 
we ar e investigating the co mpilation of command sequences into a dataflow net- 
work (jCordon aZj . l2006l) . In such a network quantum states are represented 



by a long stream of amplitudes, which has the double benefit of exposing the 
inherent parallelism while at the same time relaxing the need to fit entire vectors 
inside the same computer memory. This line of research has alrea dy lead to a 
first implementation of a parallelised simulation environment in (jVerhaegenl . 
20091) . 
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8 Appendix 

Here we give the full operational semantics of the distributed measurement 
calculus. We are concerned here with the small-step rules, indicating how atomic 
expressions in the language are evaluated. The big-step operational semantics 
of a program is found by grouping the pertaining small-step rules into one 
transition for the whole program. We use the standard notation of sequcnts 
and rules. A sequent F h _B J, is read as "given environment F, the expression 
E evaluates to the value w" . In case that the environment itself changes during 
evaluation of an expression E we write F, E — > V . Wc write F(a:;) for the value 
of X in F and F[?j/x] for the environment F with the added binding of a; to v. 
Sequents can be combined into rules which are just a different notation for 

There are four groups of rules, dealing with classical values (signals and an- 
gles) , measurement patterns and distributed measurement patterns respectively. 
Each group builds on top of the previous one. The first group of rules is to do 
with the evaluation of signals. 
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T^sum ^^^^ 

— F — ; — i ^^'^ — r — 1 — 1 TT '* 7^ j (54) 

T[v/Si\ \- Siiv T[v/Si\ h Sj I r{j) 

T \- s iv rhfj-M 
r \- s + t iv ®u ^ ' 

Here F is the outcome map or classical state, and © denotes addition in Z2. 
Angles, which can have signal dependencies percolated through via dependent 
measurements, are also purely classical. Values of signals are looked up in F via 
the ruleset for signals in order to determine the actual value of a measurement 
angle. This procedure is summarised in the following rules. 



(56) 

L r a a 

s iv t lu ^^^^ 



F h *[aY i (-!)".« + u.TT 



Fh[a]n°[a]^ ^'''^ rh*[a];*[a]0 ^^^^ 

Having defined how signals and angles are evaluated, we can now move on to 
the operational semantics of the basic commands. These commands operate on 
an environment consisting of a quantum state p as well as a classical state F. We 
have presented these rules here for density matrices; in pure state derivations we 
often use state transitions for brevity. Specifically, for a pure state we have p = 
\ip){ip\, which is mapped to L\il'){tp\L^ , with L any of the entanglement, Pauli 
or projection operators below. A pure state transition can then be alternatively 
specified as mapping to L\^p). 



(59) 
(60) 
(61) 



p,V,Eij — > AZijpAZij,T 
T \- s -Iv 

p,r,xf ^x>x-,F 

F h s ; w 
p,T,z^ z-pz-,r 



rh*[a]U/? tr(|+^)(+g|,p) 



p,r,*[Aff]-^ (+/j|.p|+/^)„r[0A] ' trp 



p,r,*[A/f]-^ (-/j|.p|-^)„r[iA] ' trp 

p,F,Ci P'X 

pS,C2C, p',F',C2 



(64) 



The first three commands are purely quantum and straightforward. The mea- 
surement command is the only command that affects the quantum state as 
well as the classical state. First, the measurement angle has to be evaluated. 
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which in turn requires evaluating the X- and Z-signals by the previous sets of 
rules. Measurement commands are also the only nondeterministic commands, 
as the measured qubit is projected onto either or \—a) with transition 

probabilities as stated. Usually, the convention is to renormalise the state after 
measurement, but we do not adhere to it here, as in this way the probability 
of reaching a given state can be read off its norm, and moreover the overall 
treatment is simpler. The last rule is for a composition of commands. 

Finally, we need a set of rules to deal with distributed extensions to pat- 
terns. Essentially these transitions describe how agents A(z,o) : Q.£ and net- 
works M = \iA.i{\i, Oi) : Qi.Si II (J evolve over different time steps. We adopt a 
shorthand notation for agents, leaving out classical inputs and output, which 
do not change with small-step reductions. 



3-1 — Aj ! Q^.Si 

a,.E ^ A, : Q,.[£,.E] 

a-" = A : Q\{q}.E ^^^^ 
a+9 = A : W {q}.£[q/x\ , 

where E is some event, and £i and £[ are event sequences. A configuration is 
given by the system state a together with a set of agent programs, and their 
states, specifically 

cr, \iTi,a.i = (T,ri,ai I r2,a2 I ... | r™,a„ . (66) 

The small-step rules for configuration transitions, denoted are specified 

below; we provide some explanations afterwards. When the system state is not 
changed in an evaluation step, we stress this by preceding a rule by a \~. 



cr,r,A : /l±)i?.[5.P] =»A cr', r',A : owi?.£: 
^2{y) = V 

o- h (ri,ai.c?a; | T2,&2-c^-V =^ ri[x i;],ai | r2,a2) 

(Th (ri,ai.qc?a: | r2, a2.qc!g =^ Ti, a+'' | r2,a2"«) 
L =^A M 



L I N =^A M I N 



(67) 

(68) 
(69) 

(70) 



Implicit in these rules is a sequential composition rule, which ensures that 
all events in an agent's event sequence are executed one after the other. The 
first rule is for local operations; we have written the full pattern instead of 
only its command sequence here to make pattern input and output explicit. 
Because a pattern's big-step semantics is given by a probabilistic transition 
system described by — >, we pick up a probability A here. Furthermore, an 
agent changes its sort depending on pattern's output O. The next rule is for 
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classical rendezvous and is straightforward. For quantum rendezvous, we need to 
substitute q for x in the event sequence of the receiving agent, and furthermore 
adapt qubit sorts. The last rule is a metarule, which is required to express 
that any of the other rules may fire in the context of a larger system. L and 
R stand for any of the possible left-, respectively right-hand sides of any of the 
previous rules, while L' is an arbitrary configuration. Note that wc might need 
to rearrange terms in the parallel composition of agents in order to be able to 
apply the context rule. This can always be done since the order of agents in 
a configuration is arbitrary. In derivations of network execution, we often do 
not explicitly write reductions as specified by (j70p . but rather specify in which 
order the other rules fire for the network at hand. It is precisely in this last 
rule that introduces nondeterminism at the network level, that is, several agent 
transitions may be possible within the context of a network at the same time. 
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